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1 Introduction 

Deep inelastic scattering (DIS) provides the primary experimental probe for determining the 
structure of the proton, which is of vital importance to the interpreations and callibration of the 
experimental program at the Large Hardron Collider (LHC). However the usefulness of DIS data 
is limited by the difficulty of a deep theoretical understanding of Quantum Chromodynamics in 
this regime, notwithstanding the progress made in the "small-x" physics from the perspective of 
perturbative QCD [1], i.e., the so-called BFKL Pomeron |2l [3l HI E] ■ The problem is that even 
in this regime the phenomenon is a subtle mixture of perturbative and non-perturbative physics. 
Recently the AdS/CFT (or more generally string/gauge duality) suggests a new approach to 
DIS at strong coupling starting from the paper by Polchinski and Strassler J6j. This article 
begins to explore the possibility of developing a new phenomenological framework based on the 
application of strong coupling Pomeron kernel computed by Brower, Polchinski, Strassler and 
Tan (BPST) 13 |9]. (For other related works, see US HI IIZI [I3-) By testing and calibrating 
the phenomenology of the BPST kernel against the DIS scattering data of HERA at small-x, 
we hope to provide a foundation for extending this approach to higher energies and to other 
diffrative processes, such as an estimation of the rate of double-diffractive Higgs production at 
the LHC. 

Before introducing the AdS/CFT strong coupling approach in Sec. pi it is useful to make some 
general observations on how we treat Pomeron physics. In AdS/CFT, non-perturbative physics 
is organized following the original observation of 't Hooft. Namely, QCD can be expanded 
(formally) term by term as a power series in l/N,. at fixed 't Hooft coupling A = gy^^c- As 
a consequence, various non-perturbative effects are classified in terms of a topological (or string 
theoretic) expansion. This has many well known qualitative successes. For example, the leading 
term for mesons is the valence approximation and for scattering Zweig-rule violating processes 
are suppressed. The nucleon is introduced as an external probe after we set Nc to its physical 
value, Nc = 3. At high energies the vacuum exchange in leading order of l/A'^c-expansion is the 
cylinder diagram, which unambiguously defines the "elementary" Pomeron as a non-perturbative 
color singlet gluonic object. This is in fact completely consistent with the weak coupling BFKL 
Pomeron, which is the leading large Nc contribution to first order in the 't Hooft coupling g'^Nc 
and all orders in ^'^A'^clog s. The BFKL equation can be viewed as the ladder approximation in 
the color singlet two-gluon channel. 

In the strong coupling limit, BPST Pomeron is now computed non-perturbatively in the 1/A^c 
expansion in leading order in A = l/gyM^c at strong coupling, with an intercept 
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2-2/VX + 0{l/X) 



Thus the two gluon weak couphng BFKL Pomeron is now viewed as Reggeized Graviton in a 
confining AdS background. In spite of the change of view point, the quahative features of weak 
and strong couphng kernels are strictly similar, due no doubt to lack of conformal breaking in 
their resepective approximations. It is instructive to compare the value of the leading power jo 
for forward scattering for both strong and weak couping as given in Fig. [T] Observe that a power 
of jo — 1.3 lies securely in the cross over region between strong and weak coupling. 
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Figure 1: The Pomeron intercept in QCD viewed as a function of the 't Hooft coupling. (A 
9ym-^c, with a = 5fy^j/47r.) This is reproduced from Fig. 12 of Ref. [7]. 



Turning to DIS scattering, (e.g. the total cross section for a virtual Compton scattering at fixed 
Q^), we now have a further tool to probe the Pomeron kernel. Varying the virtual photon's 
momentum, in BFKL language, 1/Q probes the size of the two gluon dipole, whereas from 
AdS/CFT dictionary, it probes the 5th radial coordinate, z ~ 1/Q, in the AdS space. Here we 
can ask how the BPST Pomeron can explain the observed Q^-dependence for the effective power 
behavior for virtual Compton scattering at fixed Q^ as seen in the HERA data |14j . ( See Fig. 
3 Note that the effective power increases with Q^.) 



It is important to point out that the strong coupling approach provides a natural way to include 
the non-conformal contributions due to confinement by a deformation to the AdS^ geometry. 
Thus the hard- wall AdS Pomeron provides a synthesis of the so called "soft Pomeron", i.e., a Regge 
pole which interpolates with a tensor glueball at j = 2 and the "hard Pomeron", characterized 
by the BFKL weak coupling behavior. 

The main result of this paper is to show how well the exchange of a single strong coupling 



Pomeron describes the HERA data for the small-x region of DIS. (Other recent studies include 
[l5l[l6l[I71[ig[l9l[20l[2ll|22l|23l[2l[25l[26ll2Z].) Indeed when the hard- wall model is introduced 
to implement confinement the fit is remarkably good even down to Q^ = 0.1 GeV'^. In the 
1/Nc expansion, there are of course non-linear effects which enter through eikonalization. When 
this effect becomes important for DIS, it can be interpreted as the onset of "saturation". We 
determine that for the range of energies given by HERA data these effects are still negligible for 
Q^ > 0(1) GeV'^, but that they will come into play at LHC energies. 

The organization of the paper is as follows. In Sec. |2] we begin with a short discussion on the 
kinematics for the HERA small-x regime and the standard weak-coupling partonic approach to 
DIS. We next briefly review the basic AdS/CFT formalism. In Sec.plwe treat DIS in the small-x 
limit to first order in the BPST Pomeron. We begin first with the conformal limit, and next 
discuss the modification due to confinement, using the hard-wall model as an illustration. We 
clarify how at strong coupling the small-x Regge limit and the large-Q^ limit are unified. In Sec.|4| 
we test our strong coupling results by fitting to the recently combined Hl-ZEUS small-x data 
from HERA |28) . We focus on a single- Pomeron contribution based on a local approximation 
for both the current and the proton "wave functions". We find that the confinement-improved 
treatment (hard-wall model) allows a surprisingly good fit to all HERA small-x data, with Q^ 
ranging from 0.1 GeV'^ to 400 GeV^. The single- Pomeron hard- wall fit also indicates possible 
onset of "saturation" for Q^ < 0(1) GeV"^. In Sec. ^ we carry out a nonlinear eikonal analysis. It 
is now important to fully explore the dependence of the eikonal, x(s, b, z, z'), on the 3-dimensional 
transverse space, i.e., b and z. Our analysis confirms that saturation effect is small for Q^ > 
0(1) GeV^ at HERA energy range. However, for Q < 0(1) GeV^, confinement-improved eikonal 
treatment allows a better fit to the HERA small-x data. Finally, we summarize in Sec. [6] our 
findings and discuss their possible significance. 

Note Added: After initial posting of this work, we became aware of Ref. |29) which has also con- 
firmed our prior claim, demonstrated more explicitly here, that, for the HERA range, saturation 
effect is minimal. 



2 Background 

2.1 Parton Phenomenology 

In a partonic approach, hadron structure functions Fi{x, Q"^) in DIS follow from quark and gluon 
distributions, qi{x,Q'^), qj{x,Q'^), and gk{x,Q^) respectively. At small-a;, QCD evolution die- 



tates that most of the wee partons evolve from gluons and the valence contributions become 
subdominant. At LHC, proton gluon distribution is expected to play an increasingly important 
role at small-x and/or large Q^ . In a strong coupling approach, although it is no longer mean- 
ingful to speak of quarks and gluons as the effective partonic degrees of freedom, the dominance 
of gluon dynamics at small-x justifies treating a simpler situation where quarks are first ignored, 
i.e., considering the limit N^ » Nf. In such a limit, Gauge/String duality can be applied 
directly [HC]. 




Figure 2: HERA and LHeC regions are labeled by trapezoids respectively. We have also shown 
the entire combined Hl-ZEUS small-x data points as dots. Lines for confinement and saturation 
are discussed in Sec. |3]and Sec. [5] respectively. 



In a naive parton picture, Bjorken scaling would lead to Q ^-independence for the DIS structure 
functions. In a perturbative treatment, due to QCD evolution, structure functions such as 
F2{x,Q'^) become Q^-dependent, which can be derived by using the linear DGLAP evolution 
equation. This has indeed been carried out for the HERA data, starting with phenomenological 
inputs for structure functions at low Q^, e.g., Q^ ~ 3 — 4 GeV^. Of course, an inclusive DIS 
cross section can also be treated directly as the total cross section of an off-shell photon, with 
mass squared —Q^, scattered off a proton. The Bjorken x variable is then related to the center of 
mass energy squared s of the 7*p system by x ~ ^'^/s, at high s. The limit x — )• 0, with Q^ fixed, 
can then be considered as the Regge limit for 7*p two-body scattering. The off-shell j*p total 
cross section in the limit x — )• can then be characterized by the exchange of Pomerons, as done 
for other hadronic total cross sections. This has been studied at weak coupling using the BFKL 
equation, which corresponds to summing over Reggeized gluon ladder graphs [21 |3l HI E]. It is 



therefore tempting to identify the rise of F2{x,Q'^) with 1/x, for fixed Q^, observed at HERA, 
as the onset of the Pomeron dynamics. In Fig. ^ we indicate in a log((5^) — log(l/x) plot the 
kinematic region covered by HERA and also the region of particular interest at LHC, e.g., that 
for the proposed LHeC. We also identify, by dots, the small- rr data points, with x < 0.01, for the 
combined Hl-ZEUS data set, which has been released recently [28]. There are 249 data points, 
with Q"^ ranging from 0.1 to 400 GeV'^. 

DGLAP provides a description for the evolution of the data in Q^, given Fi{x,Q*'^) at some 
initial point Q*^. The study of small-x, on the other hand, focuses on evolution in 1/x, with Q^ 
fixed. Indeed, for Q^ ranging from 0.15 GeV^ to 250 GeV^, by fitting the data to 

F2ix,Q^)^{l/xy^ff, (2) 

with Q^ fixed, "effective Pomeron" intercept, 1 -|-ee//, can be extracted jTlj . with eg// shown in 
Fig. [3] There are two intriguing aspects to this analysis. First, the effective intercept, except 
for that at small Q^, is very large, with 1 + eg// ~ 1.2 ~ 1.4. Even more puzzling is the fact 
that eeff is Q'^-dependent, contrary to a naive Regge expectation. This feature could seriously 
challenge the hypothesis of Pomeron dynamics at work in the HERA energy range. Indeed, this 
Q^-dependence has been used to support the notion of "two-Pomeron" hypothesis |30| . (For 
related earlier work, see [311 [32l [33].) 

As we have mentioned earlier, using AdS/CFT, the BPST Pomeron can provide a synthesis of 
both the "soft" and the "hard" Pomerons. Using the BPST Pomeron, we have found that the 
Q^-dependence for eg// observed at HERA, Fig. K3l can be attributed primarily to diffusion for 
Q^ large and to confinement effects for Q^ small. 

2.2 Reveiw of AdS/CFT Formalism 

In this paper, we will re-examine the DIS structure functions at small-x by treating the off-shell 
7*p total cross section using the AdS/CFT correspondence. That is, with String/Gauge Duality, 
we can provide a direct strong coupling treatment for the DIS structure functions, consistent 
with DGLAP expectation. We focus on the conformal limit as well as effects of confinement and 
evidence for the onset of saturation. We support the usefulness of this strong coupling formalism 
by testing against the recently published combined Hl-ZEUS data set. Extrapolation to the 
LHC kinematic region will also be provided. 

In our approach, two-body amplitudes are controlled by exchanging BPST Pomeron, with an 
intercept 

Jo - 2 - p (3) 
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Figure 3: Q^'dependence for effective Pomeron intercept, ap = 1 + Ce//, from [T 
curves are fits based on strong/gauge duality described in Sees. [4] and [5] 
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where we have re-written Eq. M in terms of a more convenient parameter, p = 2/vA, A = 
9ym^c being the 't Hooft couphng PJ It is useful to compare the dependence of the Pomeron 
intercept on the 't Hooft coupling at weak and strong couplings. It is worth noting that the 
BFKL Pomeron intercept for A/" = 4 SYM, calculated to second order, meets with the BPST 
intercept near ap ~ 1.2 ~ 1.3, as shown in Fig. [T] precisely in the range as extracted from 
the HERA data. At these values, the 't Hooft coupling is relatively large, around A ~ 10. 
For such a large coupling value, it is important to explore the strong coupling approach as a 
more efficient procedure for addressing physical processes such as DIS. (There have already been 
several phenomenological studies of AdS/CFT which consider A in this range. For example in 
[25], strong coupling is defined as A > 5.) 

AdS/CFT, or gauge/string duality, in principle allows a description of conformal theory at strong 
coupling by a weak coupling gravity dual in an AdS background. However, QCD can be con- 
sidered conformal only approximately at best. A conformal theory can never fully reproduce 
all experimental results due to the lack of a scale and the absence of confinement. However, 
at Q^ sufficiently large, partons inside the proton are expected to be free, and a conformally 
invariant description could be a good approximation. Conversely, at smaller Q^ values, it is rea- 



^This strong coupling result was first obtained by a direct use of AdS/CFT in Ref. [3. The fact that the 
intercept is lowered from j = 2 was strongly advocated in [34| I35j and its 1/\/A dependence was anticipated in 
[B]. This can also be arrived at by using a direct extrapolation from the weak coupling perturbative sum to the 
strong coupling limit, as done by Kotikov, Lipatov, Onishchenko, and Velizhanin [36J. See also pTj and additional 



discussion in Sec. 3.2 



sonable to expect that confinement effects should be felt. Equally important is the phenomenon 
of "saturation", which should become important due to higher order Pomeron-exchanges. In 
AdS/CFT, these non-linear effects come from eikonalization. In contrast, in weak coupling, 
saturation has been addressed primarily by considering non-linear evolutions such as the BK 
equation [38l l39l liO] . 



From AdS/CFT, high energy scattering can be visualized as taking place in a 3-dimensional 
transverse Euclidean AdS^, i.e., in addition to the usual 2-dimensional impact-space b, there is 
also the AdS-radial direction z. In [8l[9], it has been shown, for two-to-two scattering involving 
on-shell hadrons, the amplitude in an eikonal sum can be expressed as 

A{s,t) = lis f d^e'''-^ f dzdz'Pi3{z)P24{z'){l - e*x{«.M,2')}^ ^4^ 

where Pij{z) = y^—g{z){z/R)'^(j)i(z)(pj(z) involves a product of two external normalizable wave 
functions, satisfying normalization condition: f dzy^—g{z){z/R)'^(j)a{z)(t)b{z) = 6ab- The eikonal, 
X, is related to a BPST Pomeron kernel [TJ lU [9] 

x{s,b,z,z') = ^{^fK{s,b,z,z'). (5) 

2s zz 

In the conformal limit, the Pomeron kernel can be expressed in a simple closed form and the 
eikonal becomes a function of two conformal invariants, a longitudinal boost r and a variable ^ 
relating to the chordal distance in transverse AdS^, 

T = log(a's) = log{pzz's/2) , (6) 

where p = 2/vA = 2a' /B? and s = ^s. It is also possible to discuss confinement deformations, 
leading to a unified treatment for the Pomeron physics from UV to IR. For definiteness, as 
done in [71 E], effects of confinement can be illustrated analytically by adopting a model with a 
hard-wall cutoff. 

We extend in this paper the above treatment to the case involving two external currents, ap- 
propriate for DIS [6j. The structure functions in DIS can be extracted from the off-shell 7*p 
total cross sections, which, from optical theorem, are proportional to the imaginary part of the 
forward off-shell ^*p amplitudes, lra.A{s,Qi). We therefore need to consider the states of an off- 
shell photon. This will require a replacement in P13 by wave functions divergent on the AdS 
boundary. As is well known, in closed string theory there is no photon state. In [6] Polchinski 
and Strassler have considered a i?— boson propagating through the bulk that couples to leptons 



on the boundary. To be specific, for F2{x,Q), whicli corresponds to tlie sum of the transverse 
and longitudinal cross sections, the appropriate i?— boson wave functions which mimic the vector 
currents involve both Ko{Qz) and Ki{Qz). After removing the associated polarization vectors 
and one factor of fine structure constant apm , one arrives at 



1 



Pn{z) ^ Pi3{z,Q') = -{QzY{K^{Qz) + Kt{Qz)), (8) 

and, for x small where the Pomeron dynamics is applicable, 

F2{x,Q^) = ^ I dH I dz I dz'Prs{z,Q^)P24{z')Re[l-e'^^'''''''''^) (9) 

(See also |20|l41j). A similar replacement would lead to the second structure function, Fi{x, Q^), 

2xFi(x, Q^) = §^ I d'b fdzf dz'Pisiz, Q^)P2Aiz')Re (l - e^>^(^''''^'^')) (10) 

with 

Pi3{z) -^ Pisiz, Q') = ^{QzfKfiQz). (11) 

We shall next examine the usefulness of this strong coupling formalism by testing against the re- 
cently published combined Hl-ZEUS data set. To simplify the discussion, we will focus primarily 
on F2. 



3 DIS at Strong Coupling and the BPST Kernel 

Our goal is to calculate the structure function F2, which is given in terms of the cross section by 

F2{x,Q^) = -^{aT + aL), (12) 

where ax and ai correspond to the transverse and longitudinal cross sections respectively. From 
the Optical Theorem, one can express the cross section in terms of the imaginary part of the 
forward scattering amplitude, i.e., (Jtotal = {l/s)l7nA{s, t = 0). In principle Eq. Ml) could be used 
for studying any scattering process where the Pomeron exchange approximation is appropriate, 
by plugging in the appropriate wavefunctions for the states ^13(2;) and P2a{z'). For our present 
application, we consider deep inelastic scattering and the appropriate i?— boson wave functions 
for F2 can be taken to be that given by ([8|. Note that, strictly speaking, the use of ([8| is valid 
only for Q^ large and near the AdS boundary where confinement effects can be ignored. 

Treating the BPST Pomeron exchange to first order, we can then express the structure function 
F2 in terms of the eikonal linearly as 

F2(x,Q2) = (Q2/27r2) f dH f dzdz'Pr3iz,Q'^)P24iz')Imxis,b,z,z'). (13) 



with Pi^{z,Q'^) given by Eq. te|. We see that for Qz > 1, Pi3{z,Q'^) rapidly decays to zero 
because of the Bessel functions Kq and Ki. For small Qz, with ^13(2;, Q) bounded and, as we shall 
show shortly, f d?blm x{s,b,z,z') vanishing faster than 0(2;^), it follows that the integrand is 
peaked around z ~ i/q. We also note that, after a simple calculation, J„°° dz{zQ)'^Piz{z, Q"^) = 1. 

3.1 DIS in the Conformal Limit 

Let us begin by first considering the conformal limit. In this limit, the imaginary part of the 
BPST Pomeron kernel, ImK{s, b, z, z') can be expressed in a closed form, [3, and, from Eq. ([5|, 
it leads to 



lovr V vr smh ^ r / ^ 

where r and ^ are given by (pi) and M respectively. The real part. Re x{s,b,z, z'), can be 
reconstructed via a dispersion integral, or, equivalently, by solving a derivative dispersion relation. 
As noted earlier, due to conformal invariance, Imx{s,b,z,z') depends only on two variables, r, 
(or s), and ^. The 6-space integration of this kernel can be carried out explicitly [9], leading to 

2 rj cxp( ~^'°^"'~'°^^'-'^ ) 

d%Imxis,b,z,z') = f^^^(z/)e(^"^)- ^.fj -. (15) 

This expression now takes on the familiar form of "diffusion" in | log2:|. This was first noted in 
Ref. [7]. 

Let us next turn to the factor P24. For the proton, one is dealing with a normalizable state. 
Strictly speaking, this necessarily breaks the conformal limit, e.g., to model QCD we need to 
introduce a cutoff at zq ~ 1/A, where A is the confinement scale. This construct would be 
sufficient for describing glueballs. However, for baryons, one needs to go beyond the standard 
confining deformation and additional holographic prescriptions will have to be adopted |4H Ej 
I43| . It suffices to assume for now that the AdS wave- function for the proton is localized near 
the cutoff in the bulk, zq- 

Putting these two together, we have, for DIS in the strong coupling conformal limit, that the 
structure function F2 can be expressed as 

2 3/2 r l' — (log^— log^')^ \ 

F2{x, Q') = 1^ J dzdz'P,s{z, Q')P2,{z'Kzz'Q') e^^-^)- """^^ ^^fj '- (16) 



where x ~ Q'^/s and P13 given by (|8|). Note that F2 is uniquely determined once the parameter 
p and the overall strength factor Qq are specified. There are two important features for Eq. (16). 
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It is easy to verify, as noted earlier, that the integrand vanishes at least as fast as z^"^ at z = 
and is exponentially damped for z >> 1/Q. It follows that the dominant contribution to the 
integral coming from z = 0{1/Q). With z' = 0(1/Agcd), the factor e^^^^'''^ leads to 



e 



(1-P)t 



[l/xf-^ , (17) 



with Q^ fixed. The limit x — t- will lead to a fast rising in F2, violating the standard Froissart 
bound. In terms of the eikonal xi^^^j^jz'), it is generally accepted that the condition x = 
0(1) would signal the onset of "saturation". Note that this is a "local condition" in the three- 
dimensional transverse space, z and b, with z' fixed. One of the main goal of our current analysis 
is to investigate the importance of saturation in the HERA range. We will come back to this 
discussion in Sec. [H 



The second key feature for ( 16 ) is the last exponential factor 



(-{logz-logz'f\ 

exp (18 

\ pr J 

With z ~ 1/Q, this corresponds to diffusion in logQ, analogous to diffusion in "virtuality", found 
in weak coupling BFKL dynamics. As we shall see shortly, we find this diffusion effect playing a 
crucial role in understanding the Q^ dependence for the eeff{Q'^) observed at HERA, as exhibited 
in Fig. [3] 



3.2 Regge and DGLAP Connection 

Although diffusion is common to both the weak and the strong couplings, the diffusion coef- 

^2 



ficient in the strong coupling, (18), is correlated with the Regge growth, (17), leading to Q - 



independence for f^ dxF2{x,Q'^) in the large Q^ limit, i.e.. 



1 
dxF2{x,Q^) = 0. (19) 





dQ2 



It should be stressed that this feature is not shared by the conventional weak coupling BFKL 



approach. Eq. (19) is required by energy-momentum conservation, or, more technically, the 
vanishing of anomalous dimension for the second moment for the structure functions. Indeed, 
the signature feature of the BPST treatment is based on perturbing about the supergravity limit 
where the energy- momentum conservation is guaranteed. Consider the moments for the structure 
function F2, 

M„(Q2) = /" dxx"-2^2(a;, Q^) (20) 

11 



From OPE, one has, for Q large, 

M„(Q2) ^ (g)-7n (21) 

where 7^ is the anomalous dimension for the twist-two operators, appropriate for the DIS. 

From our strong coupling analysis, these anomalous dimensions can be read off the more general 
"dimension-spin" curve. Fig. [3] At strong coupling, A >> 1, one finds to leading order in 
p = 2/\/A, the relation between scaling dimension A and spin j is parabolic, 

i = 2 + p[(A- 2)2/2-1] (22) 

As explained in |7j, this follows from the physical state condition, Lq = Zq = 1, for string theory 
at strong coupling. The DGLAP dimensions and the Pomeron intercept are determined by the 
same on-shell vertex operator. The symmetry about A = 2 is due to conformal invariance, and, 
at j = 2, one has A = 4, due to energy-momentum conservation. Although we have drawn in 
Fig. [4] for weak coupling with the curve maintaining the symmetry about A = 2 and passing 
through A = 4 and j = 2, this feature is absent in a conventional perturbative treatment 



In evaluating M„, Eq. (20) acts as a J-plane Mellin transform in s, thus allowing one to identify 



the moment n with the spin j. A straight forward calculation then leads to 

7„ = 2^1 + (n - 2)/p - n (23) 

At n = 2, 72 vanishes, as expectedjj Because of our better control of the anomalous dimensions 
in the strong coupling at n = 2, our BPST Pomeron provides a dual approach to the DGLAP 
treatment, so long as a linear evolution is valid at the HERA energy. Instead of evolving in Q^, 
the x-evolution is controlled by the minimal of the A — J curve at jo = 2 — 2/\/A = 2-/9, our 
BPST intercept. 



3.3 Confinement 

Let us next turn to the case with confinement, using the hard- wall model for illustration. In 
this model, we introduce a sharp boundary in the AdS z coordinate at some value zq. This leads 
to a mass gap and a discrete spectrum of states. Our starting points are still Eqs. (12|[l3), 



but a modified Pomeron kernel has to be used. Such a kernel can formally be represented in 



^Note that, by focusing on small- 2; , our calculation for 7„ becomes unreliable for n large, which reflects strongly 
the behavior of structure functions near x — 1. 
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Figure 4: A — J curve. This figure has been reproduced from Fig. 2 of [7]. 



a spectral representation, but it cannot be expressed in a simple closed form as the conformal 



case, i.e., (14). Fortunately, for a single Pomeron contribution, we only need to evaluate the 



total contribution after integration over the impact parameter b. That is, if we go to momentum 
space, 

Xis,t,z,z')= [dHe"'^'^x{s,b,z,z'), (24) 



where t = —q'j_, we only need to know x{s, t, z, z') evaluating in momentum space at t = 0. For 
this, a closed form has been found for the BPST hard- wall kernel, |9], leading to 



Im Xhw{s, t = 0, z, z') = Im xdr, 0, z, z') + T{z, z' , r) Im Xc{t, 0, z, Zq/z'), 



(25) 



where Im Xc{t, 0, z, z') is the conformal expression given by the right hand side of Eq. ( 15 ), and 

- log ^ + pT 



T{z,z',t) = 1 — 2^fpTVTe^ erfc{r]), rj 
It follows that, with confinement, 



PT 



(26) 



327rV2 



log^ '='='/zg 



F2(x,Q^) = f^ / dzdz'Pi^{z,Q^)P2d.z'){zz'Q^) e^'-^^^ ' ^ 



1/2 



+ J-(z,z',r) 



1/2 



(27) 



The first term on the right-hand side of (27) is precisely the same as that in the conformal limit. 



( 15 ) and is model independent. The second term, which can be expressed as a linear superposition 
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Figure 5: Contour plot for coefficient function J^ as a function of log{l/ z) and log(l/x), with 



zq fixed, zq ~ Aqq^. Confinement effects become important in the region where \J^\ > 



0.3 ~ 0.5. We have also shown the entire combined Hl-ZEUS small-x data points 
by identifying Q = 1/z. 



as dots 
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of contributions from a set of image charges, is dependent upon the details of how the space is 
cutoff at zq and hence model dependent. Nevertheless it is reasonable to assume its presence can 
provide a qualitatively reliable estimate of the effects of confinement. The coefficient function, 



J-{z, z' , r), (26), at fixed z, z', goes to 1 as r — )• and to —1 as r — )• oo. Hence, at small x, i.e., 
r large, J^ — )• — 1 and confinement leads to a partial cancelation for the growth rate relative to 
that of the conformal BPST Pomeron. Moreover, since J^ is continuous, there will be a region 
over which J^ ~ 0, and, in this region, there is little difference between the hard-wall and the 
conformal results. 

In Fig. [5] we provide a contour plot for 7^ as a function of In x and In z, with z' fixed near zq. Using 
the fact that, from Pi3{z,Q'^) is peaked at 2; ~ l/Q, this can also be interpreted as a contour 
plot in 1/x and Q^, with In Q ~ — In z. Anticipating our subsequent fit which fixes the scale with 
zq ~ Aged, we have also exhibited the entire set of the recently combined ZEUS-Hl small-x data 
points as open dots on this plot. We note that, over a significant region, e.g., Q^ > 0(1) GeV'^, 
\J-'\ is small and the conformal kernel remains a reasonable approximation. Equally important is 
the observation that the region where confinement effects can be significant, roughly defined by 
the condition |J^| > c, c ^ 0.3 — 0.5, goes beyond the narrow region of Q^ = 0(1) GeV^. The 
line for c = 0.5 is also drawn in Fig. |2] indicating schematically the transition into confinement 
region. This observation goes against the conventional belief that confinement effects can be 
ignored for Q^ large, independent of x. This also casts doubt on the assumption that saturation 
at Q^ large can be understood without taking confinement into account. 



Eq. (27) is the expression which we will use shortly for comparing with experimental results. 
This expression retains the two key features noted earlier for the conformal result, namely, the 
presence of the dominant term in r, leading to (l/x)^^^ rise, and the presence of diffusion in z. 



It can also be shown that (19) still holds in the large Q limit. 



4 Confronting the HERA Data Linear Treatment 

We now carry out a test of our strong coupling DIS analysis by fitting the recently published 
combined Hl-ZEUS data set. We will first test single BPST Pomeron models, both conformal 
and with confinement, before trying the eikonal form in the next section. We restrict ourselves 
to the small-x region, i.e., x < 0.01. This remains a relatively large data set with 249 points, 
and in particular, the set extends to a large range of Q^ values, from 0.1 GeV^ to 400 GeV^. We 
find that the confinement-improved treatment allows a surprisingly good fit to all HERA small-x 
data. In contrast, the fit based on the conformal BPST Pomeron breaks down in the low-Q^ 
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region. 

Before proceeding further, we must specify more precisely the AdS wave function, (ppiz), as- 
sociated with the proton. Unfortunately, other than the expectation that it be normalizable, 
i.e., J dzy/—g{z){z/R)'^\(pp{z)\'^ = 1, one cannot determine (j)p without an explicit strong cou- 
pling model for baryon. For the current analysis, we will assume that the wave function is 
sharply peaked near the IR boundary zq, with 1/Q' < zq, with Q' of the order of the proton 
mass |4H H2| 143] . For simplicity, we will simply replace P24 by a sharp delta-function 



P24{z') « 6{z' - 1/Q'). 



(28) 



Similarly, as noted earlier, the integral over Pi3(z,Q^) in (16), and also in (27), is centered at 
z = 1/Q, with f dz{zQ)'^Pi3{z,Q'^) = 1. As emphasized earlier, (16) follows only for Q^ large. 



and there will be modifications coming from confinement for Q^ small. To simplify the discussion 
below, we will again make the local approximation by replacing Pi^{z,Q'^) by a delta-function 



Pi3iz) « C6iz - 1/Q), 



(29) 



with C ~ 1. These constitute our model-dependent parametrizations. With these substitutions, 
the z and z' integrations can be performed trivially, leading to F2 given essentially by the Pomeron 



kernel. It follows that the structure function, (16), becomes 

' / log2(Q/Q'), 

exp( ^ ^ ^ ^ 



F2{x,Q') 



9l 






pr 



(30) 



with r as a function of x and Q^: t{x, Q^) = log[(^")(;^)]- The hard-wall single-Pomeron contri- 
bution, (27), can similarly be simplified, with 77 in (26) expressed as riix, Q'^) = — ° , ^^ — -. 



To see how well our strong coupling expression comes close to the experimental data, we begin 
by first restricting ourselves to a smaller set of ZEUS data with x < 0.01 and a range of Q^ 
from 0.65 to 650GeV^. Although there are data at lower values of Q^, we do not consider them 
initially since we would like to avoid the region where confinement effect can be expected to be 
important. We will also begin by first testing against the previously published ZEUS data |14|I44). 
so we can compare the goodness of our fit relative to other published fits done by others without 
invoking AdS/CFT |45| SU U?]. Another reason for avoiding the small Q^ region is because in 
our formalism it is natural to restrict Q > Q' , where 1/Q' — z' characterizes the "size" of the 
proton and, as we shall see, will turn out to be ^ 0.5 GeV. 

We fit the F2 to the ZEUS data using Matlab with 4 free parameters, which we choose to be : p, 
go, Zq, and Q'. We have carried out fits for both the conformal and the hard- wall models. In total 
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we had 160 different data points witli 32 different values of Q^ . Surprisingly, good fits can be 
achieved for both. To be specific, we obtain the "best fit" for the confining hard-wall model, with 
the following values: p = 0.7716 ± 0.0103, g"^ = 106.01 ± 3.10, zq = 6.60 ± 1.50 GeV'^, Q' = 
0.5322 lb 0.0465 GeV . These values are "reasonable", within our general expectations. For 
instance, the value of p ~ 0.77 lies within the transition region between strong and weak coupling. 
The value of zq is also consistent with our expectation of zq ~ 0{A~^). With zq fixed, the value 
for Q' is again reasonable. The chi-square value per degree of freedom |^ we calculate for our 
best fit is x^^ r = 0.69. An equally good fit can also be obtained using the conformal Pomeron, 
with best fit values: p = 0.774 ±0.0103, g^ = 110.13 ±1.93, Q' = 0.5575 ±0.0432 GeV, and the 
corresponding chi-square value is x^o f =0-75. 

Armed with this success, we next carried out an expanded study for the recently published 
combined Hl-ZEUS data set |28| . keeping only small-x data. With x < 10~^, the set now 
extends to much smaller Q^ values, with Q^ ranging from 0.1 GeV^ to 400 GeV^, taking on 34 
different Q^ values. This is a larger data set than the ZEUS set considered above, increasing 
from 160 data points to 249 points. These are the data points shown in Fig. (Is]) as dots, with 
z = Q-K 

With the inclusion of data at Q^ < 0.65GeV^, we can no longer obtain acceptable fits by a single 
conformal BPST Pomeron. The best fit leads to an unacceptably large value of Xdo f ~ H-^- 
The fit cannot be improved by the sieve-procedure |48| . That is, one cannot attribute the poor 
fit to the presence of "outliers". In contrast, we find that the confinement-improved (hard- wall) 
treatment allows a surprisingly good fit to all HERA small-x data, especially after applying the 
sieve-procedure. We obtain the "best fit" for the confining hard-wall model, with the following 
values 

p = 0.7792 ± 0.0034, g^ = 103.14 ± 1.68, zq = 4.96 ± O.UGeV-^, Q' = 0.4333 ± 0.0243 GeV, 

(31) 
These values differ from our earlier fit only slightly and they remain reasonable, within our 
general expectations. Without sieving, we obtain X^^ r = 1.34. With the elimination of 6 data 
points as "outliers", each with x^ > 8, the X^^ c improves to be 1.16. The best cutoff is for 
X^ > 4, leading to a calculated chi-squared per degree of freedom Xdo f =1-07 for our best fit 
given above Q 

We calculated the Xd.o.f. using a standard formula, Xd.o.f. ~ jv^ Z^i -2 , where TV is the number of 

points, p is the number of parameters, and Oi,Ei and at are the observed, expected, and the total error of 
the observed values respectively. (See, for instance, P. Fornasini, "The Uncertainty in Physical Measurements", 
Springer 2008.) We did not eliminate any points in our analysis. 

*With a x^-cutoff at 8, we included a renormalization factor of 7?. = 1.0433. When the cut-off is at 4, one has 
7^= 1.2924. 
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In Fig. y we show our fits to F2{x, Q^) for all 34 different Q^ as a function of x^^. The fit for a 
single conformal BPST Pomeron is shown as blue curves and that for a single hard-wall BPST is 
shown as red curves. It is evident that, for Q^ small, the conformal model does not fit the data 
well. However, for Q^ > 4 ~ 5 GeV^, where data points exist, the conformal and the hard-wall 
models remain comparable. 

We can further test our result by turning to the HERA effective Pomeron intercept plot. Fig. 
[3] Intercept values derived from the conformal Pomeron are in red and those from the hard- 
wall model are in blue. Much of the Q^-dependence for the effective Pomeron intercept can 
be attributed to the diffusion effect. As Q^ becomes large, diffusion is enhanced, leading to an 



increase to the effect rise in F2 as 1/x increases. Let us next examine the hard- wall model, (27), 
where we see the effect of confinement is embodied in the second term. For Q^ » 1 GeV^, 
the difference between these two at HERA energy range is small. As mentioned earlier, from 
Fig. p] with |J^| << 1 mostly in this region, it is not surprising that the red and the blue lines 
converge as Q^ increases. However, as Q^ decreases, one finds \J^\ increases, and the effect of 
"destructive interference" is enhanced. In particular, in the extreme limit of z ~ z' ~ zq, the 
leading order term from the denominator of the diffusion factor, 1/\/t, largely cancel, leading 
effectively to a r"'^'^ suppression. This enhanced suppression contributes to a lowered value 
for the effective Pomeron at small Q^, which is reflected by the red-curve in Fig. [6J That is, 
using the hard-wall model as a guide, confinement leads to an effective lowering of the Pomeron 
intercept for scattering of two "soft" states, i.e., states with wave functions centered near the IR 
wall. This is consistent with the view of "soft Pomeron" dominance for hadronic cross sections in 
the near forward region, with an effective intercept ap ~ 1.08. Of course, in the region, one also 
expects eikonalization would eventually become important. We turn to this issue in the next 
Section. 

It is also worth noting that, under our local approximation where z ~ l/Q, one has, for 1/x 

large, 

r._FL^ Jo^dxx^-^\Koix)\' _ 2-p 

Ft~ f;^dxx^-p\Kiix)\'' 4-p' ^ ' 

where Fl = 2xFi and Fj- = F2 — 2xF\. With < p < 1, R ranges from 1/2 to 1/3, achieving the 
maximum in the supergravity limit when p — )• 0. This differs significantly from that expected for 
x large due to Callan-Gross relation. This is consistenlrl with that found previously in Ref. |19] . 

^In PiJ, 7? ~ (1 + Lo)/{3 + uj), with uj = 1 - p. 
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Figure 6: Global fits to the combined ZEUS-Hl sniall-x data. Dotted red lines are for single 
conformal BPST Pomeron and dotted blue lines are for single hard-wall BPST Pomeron. 
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5 Nonlinear Evolution and Saturation 



At strong coupling nonlinearity enters through eikonahzation. When eikonahzation becomes 
important, it also signifies the onset of "saturation". Instead of expanding Eq. ([9| to first order 
in the eikonal, a fully non-perturbative treatment is now required. Clearly, nonlinearity becomes 
important only when 

\x{s,b,z,z')\>0{l). (33) 

In general, both the real and the imaginary parts of the eikonal will have to be taken into account. 
Note that this is a local condition, in the three-dimensional transverse (6, z) space. (For a more 
traditional weak coupling approach, see [H |38l |39lll0] . See also |20], [26], and references therein.) 

Using the conformal eikonal as a guide, it is easy to convince ourselves that, for the HERA range, 
it is adequate to treat Pomeron linearly and one is far away from saturation. In Fig. [7^, we show 



a contour plot for the conformal eikonal, using the parameters given in Eq. (31), at zero impact 
parameter, 6 = 0. We have also made use of the fact that z ~ l/<5, which allows one to label all 
the HERA data set on the same plot. Note that, for 6 = 0, where the eikonal is the largest, the 
maximal value in the HERA range remains less than 0.2 ~ 0.3. The conformal eikonal decreases 
with increasing |6| as a power, when |6| >> z. We have carried out an eikonal analysis using the 
conformal eikonal. Not surprisingly, there is no improvement in the quality of the fit. We again 
find with unacceptable chi-square per degree of freedom, due primarily to the inadequate fit in 
the region of small Q^. In order to answer the question regarding the onset of saturation, we 
next turn to an eikonal treatment for the hard-wall model. 



I conformal at b = 



I hardwall at b = 




r, 


0.8 

y 


0.6 
1 ; ; V 


/ 


/ 








Figure 7: Contour for Im x(s, 6, z, z') for conformal and hard-wall respectively, with z' ~ zq 
fixed, as a function of x and z, (z = Q^^\ at zero impact parameter, 6 = 0. 



The eikonal for the hard-wall model, Xhw^^-, b, z, z'), can be obtained via a spectral representation 
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[U [9]. However, it is difficult to cast Xhwis, b, z, z') in a form amenable for direct manipulation, 



except for its integrated form, f d^b Im Xhw{s,b, z, z'), given by (25). Fortunately, the effect of 
eikonalization at HERA is likely to be weak except for Q^ small. We will therefore provide an 
approximate treatment which incorporates the more important new features due to confinement. 



Let us begin by first working in a momentum representation, (24) and focus on Im XhwiT, t, z, z'). 



For t 7^ 0, Im XhwiT, t, z, z') can be formally obtained by solving an integral equation 

Imxhw{T,t,z,z') = Imxhw{T,0,z,z') 

+ -4- dr' dz z'^ Imxhw{r',0,z,z)Imxhw{T -T',t,z,z') (34) 

^ Jo Jo 



with Im, XhwiT, 0, z, z') given by Eq. (25 ). Solution to this integral equation can be investigated 
numerically, which is currently under way and will be reported separately. Here an approximate 
treatment will be provided. 

The single-most important consequence of confinement is the existence of a mass gap. In the 
conformal limit, Xc^t, t, z, z') has a branch point at t = 0, which is responsible for a power-like 
fall off at large impact separation. For instance, in the supergravity limit where j — t- 2, this leads 
to the well-known cutoff for b large, 

XciT, b, z, z') ~ e^6-6 . (35) 

Incidentally, the eikonal becomes real in this limit. More generally, a power decrease signals 
the presence of a t = singularity in the momentum representation. In contrast, for hard-wall, 
because of confinement, Xc{T,t, z,z') is regular at t = 0, with its nearest singularity at t = ttiq, 
mo the mass of the lightest tensor. Because of the mass gap, confinement leads to an exponential 
damping, i.e., in the supergravity limit, it again leads to a real eikonal, but with an exponential 
cutoff 

Xhw{T,b,z,z')r^e^e-'^°' (36) 

To illustrate the effect of confinement, we show in Fig, [8] the ratio of the hard-wall eikonal to 
the conformal eikonal calculated numerically, as a function of b/zo and z/zq, with z' ~ zq fixed. 
Note the rapid drop when b > zq and, in comparison, a relative slow variation in z. 

To fully explore the consequence of confinement at finite 't Hooft coupling, it is useful to work 
with the Pomeron kernel in the J-plane. It can be shown that the eikonal for the hard-wall 
model has a cutoff at large b of the form 

Im Xhw{T, b, z, z') ~ exp[-mi& - (mo - m^f b^/Apr] (37) 
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Figure 8: Ratio of hard- wall to conformal eikonal Contour in the super-gravity limit, 
Xhwis, b, z, z')/xcis, b, z, z'). The ratio is independent of s and is plotted against z/zq and b/zo, 
with < z/zq < 1 and z' = zq fixed. 



where mi and tuq are solutions of 9^(z^ Jo(mz)) \z=zo = and dz{z'^ J2{'rnz)) \z=zo = respec- 
tively. (We have mi ~ 1.6 Zq and ttiq — 3.8 z^ . [49j). Although the J-plane is difficult 
to carry out analytically, it can be treated, for instance, numerically. For our current analy- 
sis, we shall modify our earlier the conformal hard-wall eikonal, by taking the large-b cutoff. 



(37), into account. For 6-small, as well as to take into account the proper boundary condition 
near the IR hard-wall, we shall take Im, XhwiT^b^z, z') to be of the form Im, Xhwi'^'^^'^^ ^') ~ 
Im, Xc{t^ t), z, z')+F{t, z, z')Im Xc{t, b, z, Zq/z'). We therefore adopt the following simple ansat^ 
where Im Xhwir, b, z, z') ~ D{t, b)Im Xhi(T, b, z, z'), where 



b < Zq 
b> zo 



D{T,b) — < exp[-mife-(mo-mi)2 fe2/4pT] , (38) 

Lastly, we provide an overall factor C{t, z, z') which can be fixed by the normalization condition 



(25). 



In Fig. I9I we have collected both the BPST hard- wall Pomeron fits, (dashed lines), and the 
hard- wall eikonal, (solid lines), in one place for a better visual effect. For the latter fit, real 
part of the eikonal is included dispersively. For clarity, we only show data for 17 Q^ sets, out 



This ansatz should be further improved as one moves to smaller x-value beyond the HERA range. This will 
be addressed in a future publication where we also treat elastic and total cross sections, as well as diffractive 
Higgs production at LHC. 
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Figure 9: Fit to the combined Hl-ZEUS small- a; data for F2{x,Q'^) by a hard-wall eikonal 
treatment. We have exhibited both the hard- wall single Pomeron fits, (in dashed lines), and the 
hard- wall eikonal, (in solid lines), together for a better visual comparison. The fit include 249 
data points, with x < 10~^, and 34 Q^ values, ranging from 0.1 GeV^ to 400 GeV^. Only data 
set for 17 Q^ values are shown. 
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of 34. We find that the confinement-improved eikonal treatment allows a surprisingly good fit 
to all HERA small-x data, with Q^ ranging from 0.1 GeV'^ to AOOGeV"^ and for x < 10^^, with 
a x^ = 1-04, after performing a similar sieve-procedure as done earlier. The parameters for the 
best fits remain basically the same, with 

p = 0.7833 ± 0.0035, g^ = 104.81 ± 1.41, zq = 6.04 ± 0.15 GeV~^, Q' = 0.4439 ± 0.0177 GeV. 

(39) 
Whereas these two fits are indistinguishable for most HERA data, they do begin to diverge as 
one moves to smaller x values, e.g., to the LHC range. 

Our analysis confirms that saturation effects is minimal for Q^ > 0(1) GeV^ at HERA energy 
range. For Q < 0(1) GeV^, eikonal treatment can achieve a better fit than that by a single 
hard- wall Pomeron where saturation effects can begin to be felt. In Fig. |3j we show by the 
dotted blue curve the eikonal improved effective Pomeron intercept based on hard-wall eikonal. 
Note that this further lowers the effective intercept as one moves to smaller Q^. As mentioned 
earlier, this is consistent with the view of "soft Pomeron" dominance for hadronic cross sections 
in the near forward region, with an effective intercept ap ~ 1.08. In Fig. [T]^, we provide a 
contour plots for this confinement-improved eikonal, Im Xhwis,b,z,z') at 6 = 0. By focussing 
on the eikonal at 5 = 0, we also see that it is increasing important to include non-linear effects. 



particularly for Q^ = 0(1)06^^. This of course indicates the onset of saturation, (33), and this 
is schematically represented in Fig. ^Iby a saturation line. In fact, for any Q^, for 1/x sufficiently 
small, eikonalization will always be necessary. 



6 Summary 

The AdS/CFT correspondence provides a new approach to deep inelastic scattering. In this paper 
we have carried out an analysis of the DIS structure functions at small-x using the AdS/CFT 
correspondence. Our present analysis is based on the work of Brower, Polchinski, Strassler 
and Tan (BPST) [^ where the concept of a non-perturbative Pomeron was shown to follow 
unambiguously for all gauge theories allowing String/Gauge duality. By identifying deep inelastic 
scattering with virtual photon total cross section, this allows a self consistent description at small- 
X where the dominant contribution is the vacuum exchange process. We find that the BPST 
Pomeron kernel, along with a very simple local approximation to the proton and current wave 
functions, gives a remarkably good fit not only at large Q^, dominated by conformal symmetry, 
but also at small Q^, with an IR hard- wall cut-off of the AdS. 

We first treated DIS in the small-x limit to first order in the conformal approximation limit. We 
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explain how at strong coupling the small-x Regge limit and the large-Q^ limit are unified by 
discussing the A — J curve and show how the vanishing of anomalous dimension 72 is satisfied 
automatically. We next discuss the modification due to confinement, using the hard-wall model 
as an illustration. We end with a more precise treatment of the effect of confinement on the 
structure of the eikonal, x{s, b, z, z'), in the impact space. 

This formalism is used to fit the recently combined Hl-ZEUS small-x data from HERA [28j. We 
focus on a single-Pomeron contribution based on a "local approximation" for both the current 
and the proton "wave functions". We first find that, at larger Q^, e.g., Q^ > 0(1) GeV'^, both the 
conformally invariant theory and the confined hard- wall model fit the experimental data well, e.g., 
by first restricting to a smaller set of ZEUS data |14 | I44 ] . for values of Q^ ranging from 0.65 GeV^ 
to 650 GeV^ and for x < 10^^. Armed with this initial success, we next apply our results to the 
combined Hl-ZEUS small-x data. This is a much larger data set, and, in particular, the set now 
extends to much smaller Q^ values. We find that the confinement- improved treatment (hard- wall 
model) allows a surprisingly good fit to all HERA small-x data, with Q^ ranging from 0.1 GeV^ 
to 400 GeV^ and for x < 10~^, with si- Xdo f ~ 1-07, and best fits to various parameters given 



by (31 1, e.g., with a BPST intercept at Jq ~ 1.22. In particular, we find that the Q^-dependence 
for eeff observed at HERA, Fig. K3l can be attributed primarily to diffusion for Q^ large and to 
confinement effects for Q^ small. In contrast, the conformal fits fails when the low-Q^ data is 
included. The single-Pomeron hard-wall fit also indicates possible onset of "saturation" for small 
Q^, e.g., for Q2 < o(l) GeV'^. 

Finally we carried out a nonlinear eikonal analysis. It is now important to fully explore the 
dependence of the eikonal, xi^^bjZ, z'), on the 3-dimensional transverse space, i.e., b and z. 



For the conformal limit, this is given by Eq. (14). For the hard- wall model, a more elaborate 



treatment is required. Due to a much stronger exponential cutoff in the impact parameter, 
confinement modifies drastically the conformal result. The scale of the cutoff is set by the lowest 
tensor glueball mass, which in turn depends on the confinement scale. Our analysis confirms 
that saturation effects are small for Q^ > 0(1) GeV^ at HERA energy range. However, for 
Q < 0(1) GeV^, the conformal-eikonal treatment remains inadequate. In contrast, confinement- 



improved eikonal treatment allows an improved fit to all HERA small-x data, with SiXdo f ~ ^-^^ 



and best fits to various parameters given by (39). 



Surprisingly, we find that confinement effects persist at an increasingly large value of Q^ as 1/x 
increases, as indicated schematically by a confinement line in Fig. [2] A confinement-improved 
BPST Pomeron treatment allows a surprisingly good fit to all HERA small-x data. Nonlinear 
effect due to eikonalization is small but begins to be noticeable for low-Q^ HERA data at small-x, 
indicating imminent approach of saturation. This is also represented schematically in Fig. [2] by 
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a saturation line. Note that saturation line lies above the confinement line, indicating that the 
physics of saturation should be discussed in a confining setting. Clearly, this observation is of 
significance for diffractive central production of jets, Higgs, et al. at LHC. Equally important is 
the application to the study on the onset of Froissart-like behavior |46| Wi\ [50] for DIS as well 
as the possible extension to a study on the ultra-high energy neutrino scattering and for the 
experimental search for extra-galactic neutrinos |5H [52] . Applications to these processes will be 
reported separately. Lastly, some other recent discussions related to DIS from AdS/CFT can 
also be found in[53l[ai|55l|56l|5a[58l[59l[6ni|6l]. 
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